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Since Hund's work on the ammonia molecule jl], 
the double well potential has formed a key 
paradigm in physics. Its importance is further 
underlined by the central role it plays in the Lan- 
dau theory of phase transitions [2j. Recently, the 
study of entanglement properties of many-body 
systems has added a new angle to the study of 
quantum phase transitions of discrete and con- 
tinuous degrees of freedom, i.e., spin 0, 0] and 
harmonic chains [H, @|. Here we show that con- 
trol of the radial degree of freedom of trapped ion 
chains allows for the simulation of linear and non- 
linear Klein-Gordon fields on a lattice, in which 
the parameters of the lattice, the non-linearity 
and mass can be controlled at will. The system 
may be driven through a phase transition creat- 
ing a double well potential between different con- 
figurations of the ion crystal. The dynamics of 
the system are controllable, local properties are 
measurable and tunnelling in the double well po- 
tential would be observable. 

The development of ion trap technology enables pre- 
cise control of the internal and external degrees of free- 
dom in strings of ions Q. This motivated proposals for 
the realization of massive scalar fields in ion traps [H, Q . 
These suggestions employed the longitudinal degrees of 
freedom only, leading to significant restrictions both in 
the accessible field theories, their mass etc, as well as 
the extraction of local properties of the system. We note 
that critical behavior is not experimentally accessible in 
such systems. More recently the use of radial modes in 
ion traps for the simulation of spin systems [Hi EH and 
Bose-Einstein condensation of phonons [l2| has been pro- 
posed. 

Here we suggest the use of the radial degrees of freedom 
of the ions for the creation of a scalar field whose mass we 
may adjust freely by varying the radially confining poten- 
tial. Crucially, in this way we have access to non-linear 
fields if we allow the system to approach its configura- 
tional phase transition where it changes from a linear to 
a zig-zag structure. While this transition is well known at 
the classical, high-temperature, level, HI, HI, H, 3, 12 1 
at zero temperature it realizes a quantum phase transi- 
tion in the quantum field describing the system. Further- 
more, at zero temperature the groundstate of the zig-zag 
structure is degenerate (see Fig|I]) and thus realizes a 
double well structure. Parameters such as the width and 
depth of the double well depend on the distance from the 
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FIG. 1: The geometric transition between a linear formation 
and a 2D zigzag potential. The rotational symmetry has to 
be broken by the electrode structure. By controlling the po- 
tential a double well structure could be realized for which the 
ground state is a superposition of the two zig zag configura- 
tions. 



phase transition point and may thus be controlled pre- 
cisely thanks to the high degree of control within reach in 
ion trap experiments. This approach compares favorably 
to the artificial creation of an electromagnetic double- 
well potential for a single ion by applied electric fields, 
as this would require field gradients and a degree of con- 
trol that are not likely to be accessible in the foreseeable 
future. The crucial aspect of our proposal is the fact 
that the presence of several ions in the string leads to 
the automatic creation of a well controlled double well 
potential. 

We will demonstrate that the above QPT may be char- 
acterized by the potential seen by a single normal mode, 
namely the mode of lowest frequency. Then the real- 
ization of a QPT requires the ability to experimentally 
adjust the potential perceived by the normal mode, 
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V(x) = V(0) + -a{uj)x 2 + -b(uj)x 4 + -c(w) 
2 4 o 



, (1) 



where x is the degree of freedom of the 'important' nor- 
mal mode and to is the external trap potential frequency 
analogous to the control parameter in the Landau the- 
ory. A second order transition occurs when a(u>) = 
and b(co) > 0. A First order transition occurs when the 
symmetry is broken by a cubic term. 

In the remainder of this paper we first discuss the en- 
tanglement structure of the system in both the linear 
and the non-linear regimes, i.e. next-to or away from the 
phase transition. Then we present the non-linear Hamil- 
tonian that is accessible experimentally and show how 
it enables the creation and full control of a double well 
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potential. Finally, we discuss carefully possible measure- 
ment schemes to verify our predictions. 

The Hamiltonian - In this work we consider chains of 
ions of mass m that are subject to external confining po- 
tentials in the radial (x-y) plane and axial (z) direction 
as well as their mutual electrostatic repulsion. The exter- 
nal potentials are harmonic and we will assume that the 
confining potential in the y-direction is much stronger 
than in the x-direction, resulting in an essentially two- 
dimensional problem. This may be achieved both in lin- 
ear ion traps and Penning traps (lij . Thus the Hamilto- 
nian for TV ions is of the form 
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where Xi is the distance of the i-th ion from the axis, z\ 
is the position on the axis and we have fixed the strength 
of the confining potential in the axial direction. All the 
numerical results in this work assume oj z = 1MHz. One 
may now determine the equilibrium arrangement of the 
electrostatic problem posed here. There is a critical value 
lu c (oj c w 3Nuj z /(Ay^ogN)(ior N > 1) [Hj]) of the exter- 
nal confining potentials u> x for which we observe a transi- 
tion from a linear chain (uj x > cu c ) to a zig-zag configura- 
tion (lo x < tOc) as shown in fig. [T] In the former case the 
equilibrium positions in the axial direction do not depend 
on oj x while they become functions of u> x in the latter. We 
will be interested in the regime of small displacements Xi 
when the coupling between radial and axial modes may 
be neglected [20|. Assuming the system is cooled close 
to its ground state, we may then replace the operators 
£i by their equilibrium positions which we denote by z%. 
Then we perform a fourth order Taylor expansion of the 
potential to obtain the effective Hamiltonian 
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Ki J = ~ 87 re h- Zj |s i tne coupling to axial 
modes which exists in third and fourth order has been 
omitted. While for uj x 3> lo c the contribution of the non- 
linear terms is small, this lowest order expansion of the 
Hamiltonian shows that by lowering lu x we decrease ju 
leading to an increase of (x 2 ) and thus of the contribu- 
tion of the nonlinear terms. Eventually, for lo x ~ uj c this 
will lead to a phase transition between a linear chain and 



a zig-zag configuration. Further decreasing lo x results in 
more complex spatial configurations [lit ■ 

Scalar quantum field with an adjustable mass — The 
Hamiltonian j3| describes a non-linear scalar field the- 
ory on a lattice whose effective mass may be adjusted 
via the free parameter w^Jl]]]. At the point of phase 
transition between the linear chain and the zig-zag for- 
mation (figH]) the field becomes massless corresponding 
to a critical field. The behavior of the entropy of single 
sites as well as its scaling with the size of blocks of sites 
for critical and non-critical lattice field theories has been 
of recent interest. While for non-critical systems an area 
law holds the block entanglement diverges at the 
critical point at least in Gaussian 1-D models @. Here 
we explore these questions numerically for the full non- 
linear Hamiltonian ([3]) that may be realized in ion traps 
and, in order to gain intuition, for a version of it that is 
linearized about the equilibrium positions in the linear 
and the zig-zag cases. 

We begin with the linearized model for a chain of ions, 
i.e. the regime where uj x ^> u c , and expand Xi and pi for 
the radial motion at site i into normal modes of frequency 
u! n with annihilation and creation operators a n , a^. Here 
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where b l n , are 



-*£„ bi^/hmWn/2 (a n e~ iw 
the normal mode coefficients of the crystal. Then (xf ) = 
hEn( b n) 2 h-/™w n and (p 2 ) = \ Y, n {K?timw n . At the 
phase transition between the linear and the zig-zag con- 
figuration one normal mode becomes massless, i.e. its 
frequency tu n vanishes. In a linearized model the mass- 
less mode leads to diverging (x 2 ). Since (p 2 ) remains 
finite the entropy, which is a function of (x 2 ) ■ (p 2 ), di- 
verges [1^. Only sites that have vanishing amplitude in 
the zero mode would have a finite (x ). 

For uj x < uj c , the system posseses two nearly degener- 
ate ground states and a double well potential for a large 
crystal is realized (Fig. [1]). Here the restriction to a 2- 
dimensional setting is essential as otherwise the potential 
has a 'Mexican hat' structure allowing for continuous ro- 
tations rather than forming a double well potential. For 
lj x <C <jj c we then expect that the system will possess 
slightly more than one bit of entropy per site arising from 
the binary choice of location for each ion. The excess is 
due to the entropy available to the ion within each po- 
tential well. 

At the transition point the entropy of a single ion of the 
harmonic chain diverges logarithmically. The entropy of 
a harmonic oscillator near the transition point satisfies 
Si(u x ) = log(^/(x 2 ) (p 2 )/h) « log(|6„|V(p 2 )/(2Smwo)). 
Since, the lowest normal mode frequency, to 2 = lo 2 — uj 2 , 
we find 



S 1 (u x ) « — log^/t^ - u; 2 /lu 2 z ) = 



--logu;o/w z . (4) 



For a small number of ions in the linearized model analyt- 
ical results for the entropy may be derived. For lj x > uj c 
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and two ions the entropy is S2 — — § \og(uj 2 /uj 2 — 1) + 
log ( + O (^/ujI/ujI - 1) while for the middle 
ion in a three ion chain S3 = — i log (0J 2 , /uj 2 z — W) + 
log (v/2/5/3) +1 + {{wl/ul - 12/5)V*) . 

Following the discussion of the linearized models, we 
now analyze numerically the phase transition between 
the linear chain and the zig-zag configuration for the 
non-linear Hamiltonian {3| for three ions. The transition 
between the two phases is a second order transition^ 
where the order parameter is the displacement of the 
equilibrium position from the axis as indicated by the 
Landau theory. A full numerical calculation using Hamil- 
tonian (J3j without any externally applied non-linearities 
is shown in FigO Far from the transition point u> x S> lo c 
there is excellent agreement between the full numerical 
calculation and the linearized models discussed above. 
When approaching the transition the linearized mod- 
els predict diverging single site entropies while the full 
Hamiltonian ^ yields a finite maximum value and then, 
for to x <C u> c , the entropy approaches unity i.e. an EPR 
state between one ion and the rest. The inset shows the 
comparison between the full numerics and the results of 
treating the low normal mode as decoupled, i.e., the red 
line results from the full numerics for three ions and the 
blue line is obtained by considering only the zero mode 
and then convolving the other modes to get the x\ distri- 
bution. In order to compare the results we have checked 
the equation (x?) = |((xq) + 2(x^ m ) +3(xf)), i.e., the x\ 
deviation as a function of the normal mode deviations. 
The correspondence is perfect except at the transition 
point where yj (x?) is 32nm and the decoupled normal 
mode result gives 38nm. For three ions the coefficient 
for the quartic term is b = 3 • J ■ m~ 4 which limits 
the energy gap to w 50kHz. As can be seen from Fig. 
[2] the nonlinear coupling between the different normal 
modes is negligible. As the number of ions increases the 
distances between the ions decreases; this raises the min- 
imal energy gap and reduces the maximal entanglement. 
In this region the chain becomes highly non linear; this 
enables the observation of features peculiar to nonlinear 
systems like solitons and quantum and classical discrete 
breathers [2^]. 

After this discussion of the critical entanglement prop- 
erties at the phase transition between the linear chain 
and the zig-zag formation we now move to the discussion 
of the second main feature of this system. 

Double well potential - The formation of double well 
potentials in ion traps through the application of external 
potentials has been considered in the past as an interest- 
ing system for the exploration of the classical-quantum 
boundary and as a system in which to preform atom in- 
tcrferomctry. A number of groups are currently working 
towards scaling up ion trap quantum information pro- 
cessing to larger numbers of qubits following the general 




1.2 1.3 1.4 1.5 1.6 1.7 1.8 



FIG. 2: Transition of three ions. The green curve shows the 
entropy of the first ion as a function of uj x . At the transition 
point Uc/wz = (12/5) 1//2 the entropy in the Gaussian approx- 
imation diverges and the real value (with no approximations) 
reaches a maximum value after the transition point. The in- 
sets show the probability distribution before (right) and after 
(Left) the transition. The full numerics is in red and the re- 
sults obtained from considering only the normal mode are in 
blue. The distance is in 4.5nm units. 



approach suggested by [28(. A key step in this approach 
is to separate two ions that are initially close together in 
a particular trapping zone of an array of miniature traps, 
so that they end up in separate trapping zones. This is 
done by raising the potential on a small 'separation' elec- 
trode. This electrode needs to be very small in order to 
allow for the production of a steep gradient in the po- 
tential between the ions. During the separation protocol 
the trap transiently moves through a double well config- 
uration. The separation of a pair of ions in this way has 
been reported using a separation electrode with an ax- 
ial dimension of lOO^m [29J. A new generation of traps 
with smaller dimensions is currently under development 
worldwide but problems with anomalously high heating 
rates in these very small traps have yet to be solved. A 
theoretical study of relevance to this problem indicates 
that higher field gradients may be generated with rela- 
tively large electrodes that provide a dc octupole moment 
along with the usual a.c. radial quadrupole [3fj| . Traps 
based on these ideas have yet to be tested. Despite these 
advances in the experimental ability to realize double well 
potentials for ions with a measurable tunneling rate is 
still not within reach. However, as we have explained 
above, the zig-zag configuration gives rise to degenerate 
ground states separated by a tunnel barrier, thus provid- 
ing a natural double well potential requiring no externally 
applied fields. As only the harmonic trapping potential 
needs to be controlled, precise control of the depth and 
separation of the effective double well potential and mea- 
surable tunneling rates should be achievable. This moves 
the realization of double well experiments in ion crystals 
within the reach of current experimental technology 
We will now discuss various possibilities to control the 



parameters of the double well, to rotate the two-level 
system formed by the double well potential across the 
entire Bloch sphere and to measure coherence and tun- 
neling rates in the double well potential. In the vicin- 
ity of the transition point the low normal mode alone 
controls the transition. Therefore by controlling the pa- 
rameters of the Hamiltonian ([3]) and thus the values of 
a and b in EqfTJ a double well structure with a distance 
between the two minima is 2(-a/b) 1/2 and the width of 
the wave function of (H/ (2^— ma)) 1 / 2 is created. The 
double well is created at the stage when each well is deep 
enough to accommodate one level, i.e., \\ ~ h^/2a(the 
depth is of the order of one excitation). Which means 
a > 2 5 / 3 (h 2 b 2 1 m) 1 / 3 . For the zig-zag configuration for 
three to ten ions y^ui 2 . — lo 2 varies between 76kHz to 
155kHz and for the square configuration it is 36kHz. 

Such double wells can be created for the zig-zag con- 
figurations. In order to make the transition between the 
single well and the double well adiabatic, the transition 
rate should be slower than the minimal energy gap. The 
energy gap that should be considered is of the order of the 
energy gap of the quartic potential. The ground state in 
the Gaussian approximation is ^gyst^^b 1 / 3 and the en- 
ergy gap is approximately (|) 4 3 Yfil^b 1 ! 3 . In the zig zag 
configuration the minimal energy gap varies from 80kHz 
for three ions to 200kHz for seven ions. High energy gaps 
increase the robustness of the system to decoherence at 
the transition point. 

Double wells can also be realized in the transition be- 
tween 2D and 3D. This transition may be more suitable 
for Penning trap crystals which arc not heated by micro- 
motion. The zig zag transition could be realized both in 
Paul traps and in Penning traps after axialization. In the 
case of four ions, for uj x /lu z < 0.822 the stable configura- 
tion is two dimensional in the x — y plane (FigJ3Ja)). For 
0.822 < lo x /lo z < 1.27 a tetrahedron shown in Figj3^b) 
is the stable configuration. In the square configuration 
b = 5 ■ 10~ 4 J ■ to thus the energy gap is 36kHz and 
the local potential is 140kHz at the optimal double well 
point. 

The tunneling rates and the minimal energy gap in- 
crease with the number of ions. For three ions in the zig 
zag configuration the optimal point is at ui x = uj c — 3kHz 1 
where distance between the configurations is 180nm and 
the tunneling rate is 26.1kHz. For seven ions the opti- 
mal point is at u> x = u> c — 5.1kHz the distance between 
the configurations is ~ 1/im and the tunneling rate is 
55kHz. For the square configuration the optimal point 
is at oj x = u> z — 1kHz and the distance is ~ 300nm. For 
the zig zag configuration the tunneling rates vary from 
3kHz for three ions to 4.5kHz for seven ions and 1kHz 
for the square configuration FigJSJa). The distance at the 
optimal point varies between 30nm to AOnm and 20nm 
for four ions. 

Having created a double well potential in the way 
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FIG. 3: (a) ID Square configuration of ions, (b) is a tetra- 
hedron phase, (c) From single well to double well adiabatic 
passage. Final fidelity as a function of rate. In order to correct 
this error, cooling in the double well phase could be possible. 
In this phase the system would choose randomly the left or 
the right state, (d) Rabi flipping for four ions 



described above the manipulation of its quantum state 
may be achieved in two ways. The first approach maps 
the external degree of freedom into the electronic one 
employing a Hamiltonian r)Q.i(af om \L - R) (L + R\ + 



a at om \l + R) (L - R\). This is achieved employing a 
laser driving an electronic transition with a detuning 
equalling the energy gap between \L — R) and \L + R) 
and is possible since in the regime where each well sup- 
ports one localized state and the spatial separation is less 
than a [im. one laser can make the transitions as long 
as the operation rate is slower than the energy gap figH] 
(10kHz for three ions). Thus the manipulation of the ex- 
ternal degree of freedom is achieved by mapping it onto 
the electronic degree of freedom, followed by its manip- 
ulation employing an on-resonance laser and the subse- 
quent mapping back onto the external degree of freedom. 

For larger separation two lasers could create Raman 
transitions between the \L) and |i?)(the ground states of 
the the wells) states using an excited state \C) of the 
system whose wave function extends across both wells 
and whose energy is thus well above the gap between 
the two lowest lying states of the double well potential. 
The Raman transition could be created using the inter- 
nal degrees of freedom of the ion, by using the three state 
|| L) < — > |t C) < — ► | J. R). Raman transitions for dou- 
ble well systems were suggested in [31, 32, 33 1. For this 
purpose an asymmetry between the two potential wells 
should be introduced in order to distinguish between the 
right and the left well which will break the symmetry be- 
tween left and right. In order for the cubic term(ax 3 ) to 
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FIG. 4: The Double well energy levels for the three ions case. 




FIG. 5: Cut-away view of a cylindrical Penning trap con- 
sisting of 5 electrodes showing equipotentials calculated using 
SIMION. With the electrodes labelled a-e from left to right 
the applied voltages are: V a = 15.8V, V b = -3V, V c = OV, 
Vd = 3V and Vd = 4.5V . The inner diameter of the electrodes 
is 6mm. The widths of the three inner electrodes (in the axial 
direction) are 0.8mm and the voltages applied to them have 
been chosen to generate a significant axial cubic term. The 
voltages on the outer electrodes are chosen to push the po- 
tential minimum back to the middle of the structure and to 
provide the usual quadratic term. The inset shows the axial 
potential in the central region of the trap. 

create an energy gap of 1kHz, a should be of the order 
of lO -10 -^-. This will require application of 10 V voltage 
for electrode size of 1mm. A cubic term can be created 
using a trap geometry shown in Fig. [5j 

The manipulation of the external degree of freedom 
may also be achieved by applying a radio frequency drive, 
i.e., by applying a time varying potential Vq cosut. This 
will create the effective a x or a y rotations. Another way 
to create only one of these rotations is by adding a small 
cubic interaction which will add a phase shift between 
the \L) and \R) states. 

An alternative approach to the manipulation of the 
state of the double well system exploits the fact that 
we are able to control the state of a single well poten- 
tial. Thus by moving back and forth across the phase 



transition between single well and double well potential 
adiabatically we can generate any desired state in the 
double well potential by controlling the state in the sin- 
gle well potential. The state |0) of the single well would 
change adiabatically to \L) — \R) in the double well while 
the state |1) evolves into \L) + \R). Due to the non- 
linear term - bx 4 , the minimal gap would not vanish 
in this process and thus the transition can be carried 
out at a sufficiently high rate while remaining approxi- 
mately adiabatic. Thus, in order to move adiabatically 
from one phase to the other, the rate should be less than 
the minimal energy gap but higher than the decoherence 
rate (FigGHc)). Arbitrary measurements on the double 
well may be achieved by combining unitary rotations im- 
plemented as explained above with measurements of the 
presence or absence of the system in the left or right well. 

In the following we show a specific example of the tran- 
sition between the two phases. For four ions, in the tran- 
sition between a 2D and a 3D structure, the adiabatic 
condition is restrictive at the transition point luq = 0, at 
which the potential is V = jbx 4 . For this transition the 
the minimal energy gap between the ground state and 
the first excited state is 50kHz. This implies that the 
adiabatic sweeping could be relatively fast and create a 
small number of excitations. Fig. [3] shows the overlap of 
the state at a specific time of the transition, with the in- 
stantaneous ground state. The faster the rate the smaller 
the overlap. Fig. EJd) shows the Rabi flipping for four 
ions between the different states in Fig. OJb). It can be 
seen that a few cycles could be measured. Another way 
to measure the state is to change the system adiabatically 
to the state where the distance between the two wells is 
larger than one wavelength and then to measure, the \L) 
and \R) states using lasers. 

The realization of the double well described here could 
open the way to interferometry at the Nano-scale. Ion 
traps have already been used as a measurement appara- 
tus at the nanoscale [IH ■ Double wells for single atoms 
and BEC's have been used for precise measurement in 
interferometry experiments. Using the Ion trap double 
well, measurements of non linear electric fields on the 
nanoscale should be possible. Measurements via inter- 
ferometry should be more precise than the single well 
measurements. Moreover since linear fields are not seen 
by the ion trap, this procedure could measure cubic 
fields while completely eliminating the linear contribu- 
tion. This double well could also measure magnetic field 
gradients. The initial state of the internal degrees of 
freedom can be initialized in a superposition of two lev- 
els and the phase induced by the magnetic gradients can 
be measured by the population of the excited state. 

The ideas presented here may be realized in linear Paul 
traps as well as Penning traps. Linear ion traps will 
exhibit micro-motion when ions are displaced from the 
trap axis. Nevertheless, for small deviations as discussed 
here, micromotion and the resulting decoherence effects 
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will not be significant. The advantage of the Penning trap 
is that it does not suffer from micro-motion, and thus 
the crystal may be less fragile to decoherence. However, 
observations need to take place in a rotating frame in 
which the effect of the magnetic field is cancelled [24j ]. 

DISCUSSION 

We have suggested a realization of a non-linear Klein- 
Gordon field on a lattice with adjustable parameters. For 
a small number of ions the crystal may be cooled to the 
ground state and entanglement properties may be mea- 
sured to a high precision. We have demonstrated that by 
crossing the quantum phase transition between a linear 
and a zig-zag configuration a fully controllable double 
well potential may be achieved with available technol- 
ogy. This demonstrates the potential of this system for 
the study of complex quantum many body phenomena. 
By controlling the nonlinear terms for example, highly 
non - linear excitations, such as kinks and breathers, 
can be created. Furthermore, rapidly driving the sys- 
tem through the phase transition, the Kibble mechanism 
may be observed [36| , due to the ability of of addressing 
and measuring individual ions with high precision. 
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